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ABSTRACT 

We have conducted a spectroscopic survey of the inner regions of the Sagittarius (Sgr) dwarf galaxy 
using the AAOmega spectrograph on the Anglo-Australian Telescope. We determine radial velocities 
for over 1800 Sgr star members in 6 fields that cover an area 18.84 deg 2 , with a typical accuracy of 
w 2kms~ 1 . Motivated by recent numerical models of the Sgr tidal stream that predict a substantial 
amount of rotation in the dwarf remnant core, we compare the kinematic data against N-body models 
that simulate the stream progenitor as (i) a pressure-supported, mass-follows-light system, and (ii) a 
late-type, rotating disc galaxy embedded in an extended dark matter halo. We find that the models 
with little, or no intrinsic rotation clearly yield a better match to the mean line-of-sight velocity in 
all surveyed fields, but fail to reproduce the shape of the line-of-sight velocity distribution. This 
result rules out models wherein the prominent bifurcation observed in the leading tail of the Sgr 
stream was caused by a transfer from intrinsic angular momentum from the progenitor satellite into 
the tidal stream. It also implies that the trajectory of the young tidal tails has not been affected by 
internal rotation in the progenitor system. Our finding indicates that new, more elaborate dynamical 
models, in which the dark and luminous components are treated independently, are necessary for 
simultaneously reproducing both the internal kinematics of the Sgr dwarf and the available data for 
the associated tidal stream. 
Subject headings: 



1. INTRODUCTION 

Sercndipitously discovered in a kinematic survey of the 
outer regions of the Galactic bulge, the Sagittarius (Sgr) 
dwarf galaxy (Ibata et al. 1994) displays one of the clear- 
est cases of hierarchical galaxy formation in action, as it 
is currently being cannibalized by the Milky Way's tidal 
field and is close to full tidal disruption (Velazquez & 
White 1995; Ibata et al.1997; Niederste-Ostholt et al. 
2010). The existence of an associated stellar stream that 
fully wraps around the Galaxy (Mateo et al. 1996; Ibata 
et al. 2001; Martmez-Delgado et al. 2001, 2004; Dohm- 
Palmer et al. 2001; Majewski et al. 2003; Belokurov 
et al. 2006) provides clear-cut evidence that the tidal 
stripping of stars commenced several orbital periods ago 
(Edelsohn & Elmegreen 1997; Helmi & White 2001; Law 
et al. 2005; Fellhauer et al. 2006). 

However, in spite of the extensive theoretical efforts 
devoted to reproducing the position and velocity of the 
Sgr system (remnant core plus stellar stream) , there cur- 
rently is no model able to describe the wealth of existing 
observational data. A number of aspects have proven 
particularly difficult to understand. First, the width and 
velocity dispersion of the detected parts of the stream 
suggest that the stellar material was tidally stripped in 
a relatively recent epoch (~ 2-3 Gyr ago). This seems 
at odds with the modelled orbit of the progenitor dwarf, 
which has an estimated perigalacticon of ~ 15 kpc and 
a period ~ 1 Gyr (Ibata et al. 1997; Zhao 1998). Jiang 
& Binney (2000) showed that the relatively young age of 
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the stream could be explained if the progenitor mass was 
originally as large as ~ 10 11 M Q and its accretion onto 
the Milky Way (MW) occurred at a large Galactocentric 
distance, jsj 200 kpc. In this scenario, dynamical friction 
would be responsible for dragging the Sgr dwarf to its 
current orbit on a time-scale that spans several orbital 
periods. 

Alternatively, motivated by the discovery in Sloan Dig- 
ital Sky Survey data (SDSS; York et al. 2000) of a bifur- 
cation in the leading tail of the Sgr stream (Belokurov 
et al. 2006), Fellhauer et al. (2006) explored numeri- 
cal models where the original progenitor mass was small, 
10 8 -10 9 M Q . In these models the role played by dynam- 
ical friction becomes inconsequential. As a result of the 
null orbital decay, the mass-loss process spans over a 
large number of orbital periods, which leads to the for- 
mation of multiple stream wraps. If the shape of the 
MW Galactic potential is close to spherical, these wraps 
may appear as bifurcated streams on the sky. Recent 
estimates of the stellar mass in the known pieces of the 
Sgr stream have lent some support to this scenario, as 
the coadjuted brightness of the Sgr system seems com- 
parable to those of the M31 dwarf ellipticals (dEs) NGC 
147 and NGC 185 (Niederste-Ostholt et al. 2010), whose 
dynamical masses are <~ 1O 9 M (Geha et al. 2010). 

However, detailed analysis of SDSS observations sug- 
gest that both arms of the bifurcated stream may exhibit 
similar distances, velocities and metallicity distributions 
(Yanny et al. 2009; Niederste-Ostholt et al. 2010), which 
is hard to accommodate by models that explain the bifu- 
raction as two wraps of different ages, or as independent 
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Fig. 1. — Upper panel: the spatial distribution of selected 2MASS 
candidate Sgr K and M-giants (see text for details). In red (bold 
points) we highlight stars with measured spectra from 6 fields (Fl, 
F6) of 2 degree diameter. The large concentration of stars to- 
ward the upper right denotes the main Milky Way contamination. 
Lower panel: the GSR velocity distribution of the candidates ob- 
served. The dashed line at 120 km s -1 is our selection boundary 
for Sgr members and provides a clean separation from Milky Way 
foreground interlopers. 



streams from different progenitors. 

In a recent contribution Penarrubia et al. (2010, here- 
after P10) were able reproduce the existing data on the 
leading tail bifurcation by adopting N-body models in 
which the Sgr dwarf galaxy was originally a late-type, 
rotating disc galaxy embedded in an extended dark mat- 
ter halo, rather than a non-rotating, pressure-supported 
dwarf spheroidal galaxy (dSph), as previously thought. 
In this scenario, the two arms are both leading the galaxy, 
and arise naturally during consecutive pericentric pas- 
sages, which occur at different orientations of the Sgr disc 
relative to its orbit. Furthermore, these authors showed 
that internal rotation can efficiently alter the trajectory 
of the stream with respect to the orbit of Sgr, which 
might indirectly affect any constraint on the shape of the 
MW halo derived from pressure-supported dwarf models. 
It is thus clear that understanding the bifurcation has im- 
portant implications for the Galactic potential, which has 
been inferred from the Sgr stream to be variously spher- 
ical (Ibata et al. 2001), oblate (Johnston et al. 2005), 
prolate (Helmi 2004), or triaxial (Law & Majewski 2010, 
hereafter LM10). 

Clues to the true origin of the Sgr dwarf (dSph, dE or 
late- type spiral galaxy), as well as to the shape of the 
MW dark matter halo, can be gleaned from the internal 
kinematics of the remnant core. Here we present the 
results of a spectroscopic survey of the Sgr remnant core 
using the AAOmega instrument at the Anglo- Australian 
Telescope (AAT). 

2. OBSERVATIONS 

We observed candidate Sgr red giant stars with the 
AAOmega fiber-fed spectrograph on the AAT the nights 
of 2010 June 3-6. AAOmega obtains spectra of ~ 400 
targets simultaneously over a 2-degree field of view, em- 



ploying a dichroic to send light to separate blue and red 
spectrographs. We used 1500V (R ~ 3700) and 1700D 
(R ~ 10000) gratings in the blue and red, to cover the 
regions around Mg I 5170A and the Ca II 8500A triplet, 
respectively. Although most of the observing run was 
lost to clouds, we were able to obtain high-quality data 
on 6 out of 9 intended fields in the central regions of the 
Sgr dwarf. The data were reduced using the Australian 
Astronomical Observatory's 2df dr software package. 

Target stars were selected from the 2MASS PSC cat- 
alogue (Skrutskie et al. 2006). The photometry was 
first reddening-corrected using the Schlegel et al. (1998) 
extinction maps, and then candidates were required to 
satisfy the criteria 9 < K B < 13 and J ~ K B > 2.222 - 

0. 111-Ks, i.e., likely Sgr K and M-giants. The location 
on the sky of these candidate Sgr stars is shown in the 
upper panel of Fig. 1; objects for which we obtained 
spectra are shown using bold points. The lower panel 
of Fig. 1 shows the measured Galactic Standard of Rest 
(GSR) velocity distribution as a function of RA. Veloci- 
ties were derived from cross-correlating the Ca II triplet 
region of the spectrum (8450A- 8750A) with a model 
template spectrum (see Battalgia et al. 2008 for further 
details), and were corrected to a GSR frame. Star-by-star 
comparison of these results with an independent cross- 
correlation of the spectra, using the IRAF task f xcor 
and a synthetic Ca II triplet spectrum, yielded mean off- 
sets of J5 0.5 km s" 1 , with a typical standard deviation of 
< 1.5 kms -1 . 

MW foreground contamination is readily recognized 
from the spatial density variations visible in the upper 
panel coupled with the clean separation in a GSR frame 
from Sgr members (wgsr > 120 kms" 1 ) that appears in 
the lower panel. We find that 1805 stars from our spec- 
troscopic sample satisfy the velocity cut and are consid- 
ered as Sgr dwarf members in this paper. 

3. THEORETICAL MODELS 

To analyze whether the Sgr dwarf possesses internal 
angular momentum we compare the above kinematic 
data against two types of N-body models: one in which 
the Sgr dwarf galaxy is composed of a late-type expo- 
nential disc embedded in a non-singular isothermal dark 
matter halo (rotating model), and another where the 
stellar and the dark matter particles follow a spherical 
King (1966) model (pressure-supported model). The main 
characteristics of the rotating model are the following (for 
details see P10): the initial mass is M = + — 
3.5 x 10 8 M Q + 2.4 x 10 9 M Q » 2.8 x 10 9 M Q , where m d 
and nih denote the disc and dark matter halo masses. 
The initial disc scale-length and the halo core radius are 
~ 0.9 kpc and r^ = 0.5Rd, respectively. The rota- 
tion curve of this model peaks at « 43 kms" 1 . For the 
pressure-supported model we adopt similar parameters 
as those derived by LM10 in order to reproduce most of 
the existing observational constraints on the Sgr stream, 

1. e. an initial mass of Mq = 6.0 x 10 8 M Q , a core radius 
Rk — 0.85 kpc, and a King tidal radius R t = 5.02 kpc. 

We assume a Galactic potential composed of disc, 
bulge and dark matter halo (see P10 for details). A 
brief description of this potential is the following: the 
disc component follows a Miyamoto-Nagai (1975) model 
with a mass Md — 7.5 x 1O 1O M0, and radial and ver- 
tical scale lengths a — 3.5 kpc and b = 0.3 kpc. The 
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MW bulge follows a Hernquist (1990) profile with a mass 
M b = 1.3 x 10 10 M Q and a scale radius c = 1.2 kpc. The 
MW dark matter halo is modelled as a Navarro, Frenk & 
White (1997) profile with a virial mass M vir = 10 12 M o , 
virial radius r V i r = 258 kpc and concentration c V i r = 12 
(Klypin et al. 2002). In order to construct a Sgr 
dwarf model that matches the position and velocities of 
the observed stream pieces, we follow LM10 in assum- 
ing that the Galactic halo is triaxial in shape. This is 
done by introducing elliptical coordinates with the sub- 
stitution r — > m, where m 2 = x 2 /a 2 + y 2 /b 2 + z 2 /c 2 , 
and (a, 6, c) are dimensionless quantities with values 
(a, b, c) ~ (0.61,1.34, 1.22). 

We use test particles to integrate the orbit of the 
Sgr dwarf back in time from the current position and 
velocity of the remnant core in order to derive initial 
conditions for our N-body models. For simplicity, we 
hold the halo parameters fixed through the evolution 
of our Sgr models, given that tidal streams are barely 
sensitive to the past evolution of the host potential 
(Penarrubia et al. 2006). We use an total integration 
time of 2.5 Gyr (approximately 2.5 orbital periods), as 
this is the (minimum) time required to reproduce the 
observed pieces of the Sgr stream (see P10). The cur- 
rent position of the Sgr dwarf in galactocentric coor- 
dinates is (£>,/, 6) S gr = (25kpc,5°.6,-14 .2). The ve- 
locity vector can be derived from the line-of-sight ve- 
locity vios — 137 km s _1 (Ibata et al. 1997), and the 
proper motion measurements of Dinescu et al. (2005), 
who find cos b,/i b ) = (-2.35 ± 0.20,-2.07 ± 0.20) 
mas/yr. Adopting standard values for the Galactocen- 
tric position of the sun (X,Y,W) Q — (—8,0,0) kpc, 
the Local Standard of Rest and the solar peculiar mo- 
tion (U,V,W) & = (10, 225, 7) km s" 1 (Binney & Mer- 
rifield 1998), we find a Galactocentric space motion is 
(II, 6, WOsgr ^ (HO, -37, 264) km s~ l . After integrating 
our N-body representations of the Sgr stream progenitor 
forwards in time for 2.5 Gyr, the rotating and pressure- 
supported models have respectively lost 42% and 48% of 
their initial stellar masses to tides. 

Fig. 2 shows the stellar remnants of our rotating (up- 
per panel) and pressure-supported (lower panel) dwarf 
galaxy models. N-body particles are color-coded accord- 
ing to their mean line-of-sight velocity in the GSR frame. 
This figure illustrates a few interesting points. First, 
both dwarf models show velocity gradients throughout 
the dwarf remnant core. Indeed, this is a feature ex- 
pected in models where the Sgr dwarf is assumed to 
have originally been a late-type galaxy. In particular, 
the upper panel of Fig. 2 shows that, although tidal mass 
stripping removes a large fraction of the original angu- 
lar momentum in the progenitor Sgr disc, the remnant 
core still rotates with a velocity amplitude ~ 20 km s _1 , 
which translates into a net velocity difference of ~ 40 
km s _1 between stars at declinations above and below 
5~ -30°. 

The presence of a velocity gradient in the pressure- 
supported model, however, may be a result of the line- 
of-sight projection of the orbital velocity vector, rather 
than a manifestation of internal rotation. Since the or- 
bital trajectory (dashed line) is practically perpendicular 
to the MW disc (dotted line), we expect a velocity gra- 
dient to principally manifest at different latitudinal di- 
rections with an amplitude Av app = Ab(dr/db ■ v) ~ 
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Fig. 2. — The remnant cores of the Sgr dwarf rotating (upper 
panel) and pressure-supported (lower panel) models. Particles are 
color-coded according to the mean line-of-sight velocity. Circles 
denote the location of the 2dF camera fields. Dashed lines and 
arrows mark the orbital trajectory of the Sgr dwarf. The location 
of the Galactic equator is shown with a dotted line for ease of 
reference. Note that the remnant core of a pressure-supported 
model shows a velocity gradient with an amplitude ~ 10 km/s 



— 231A6kms~ 1 , where Ab = b — &g gr . Given that 
the maximum separation of our fields (white circles) is 
| Aft | ~ 4°, we expect the line-of-sight velocity to in- 
crease/decrease approximately 16 kms -1 with respect 
to the systemic velocity of the dwarf in the fields that lie 
closest to/furthest from the MW plane. 

The second interesting point refers to the direction of 
the velocity gradient. Whereas the rotating model pre- 
dicts an increasing velocity towards higher Right Ascen- 
sions, the pressure-supported model shows the opposite 
trend. Given the accuracy of our velocity measurements 
(typically 2 kms -1 ) and the large number of Sgr mem- 
bers in our sample (N w 1800), the above model predic- 
tions are straightforward to test. 

4. MODELS VS. DATA 

Fig. [3] shows the comparison of the velocity distribu- 
tion expected for rotating (dashed lines) and pressure- 
supported (dotted lines) dwarf models in each of the 
fields where data (solid lines) exists. A number of re- 
vealing points can be gleaned from this figure. The 
first is that we observe a clear variation in the veloc- 
ity peak along the major axis of the dwarf. For exam- 
ple, comparing fields Fl and F4, which are separated 
by 8 degrees (see Fig. 1), the velocity distribution peaks 
at ~ 185 km s -1 and ~ 165 km s -1 , respectively. The 
amplitude, as well as the sense, of the velocity shift is 
consistent with the variation in the mean line-of-sight 
velocity expected from a projection effect, as previously 
discussed. Further evidence for this interpretation can be 
found in the absence of a mean velocity gradient along 
the minor axis, as the peak velocity in both the F2 and 
F6 fields is 175 km s -1 . 

Indeed, it is clear that the pressure-supported model 
provides a better match to the mean velocities mea- 
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Fig. 3. — Line-of-sight velocity (in Galactic Standard of Rest 
coordinates, see text) histograms in our 6 AAO fields (solid lines) in 
bins of Av = 4 km/s (note that the Sgr member stars have velocity 
measurements with uncertainties e v < Av). For ease of comparison 
we show the predictions from pressure-supported (dotted lines) and 
rotating (dashed lines) models. Note that the pressure-supported 
models provide a better fit to the data. 



sured throughout the Sgr remnant core than the rotating 
model, which predicts a velocity trend that is nearly op- 
posite to the one observed. Our data thus seem to rule 
out the scenario proposed by P10, wherein the bifurca- 
tion observed in the Sgr stream would be caused by a 
transfer of internal angular momentum from the progen- 
itor satellite into the stream. 

However, it is worth noting that, although the 
pressure-supported model provides a reasonable qualita- 
tive description of the mean velocities observed through- 
out the remnant core, the shape of the line-of-sight ve- 
locity distribution appears considerably more leptokurtic 
than predicted in all the surveyed fields, except perhaps 
in Fl. This result is hardly surprising, considering that 
the dwarf model adopted by LM10 in order to reproduce 
the properties of the Sgr stream assumes for simplicity 
that the stellar and the dark matter particles follow the 
same density profile. Our data clearly indicate that bet- 
ter kinematic coverage of the Sgr core plus more sophis- 
ticated mass modelling should provide strong constraints 
on the original distribution of stars and dark matter in 
this galaxy. 

5. DISCUSSION 

Despite the prodigious theoretical efforts undertaken 
to reproduce the characteristics of the Sagittarius tidal 
stream (e.g. Ibata et al. 2001; Martfnez-Delgado et al. 
2004; Helmi 2004; Law et al. 2005, 2010; Fellhauer et 
al. 2006), two aspects of the stream have proven partic- 
ularly challenging to understand. First, the leading tail 
of the Sgr stream is bifurcated, with both arms exhibit- 
ing similar distances, velocities and metallicity distribu- 
tions (Yanny et al. 2009; Niederste-Ostholt et al. 2010), 
which would appear to refute earlier models that explain 
the bifurcation as two wraps of different ages (Fellhauer 
et al. 2006), or as independent streams from different 



progenitors. Second, the position on the sky and the he- 
liocentric velocities of stream members suggest the possi- 
bility of the MW dark matter halo being triaxial in shape 
(LM10). However, the orientation of the favored triaxial 
halo model is puzzling, for its intermediate axis appears 
to be aligned with the spin vector of the MW disc. This 
is hard to understand, as circular orbits about the in- 
termediate axis are unstable, raising questions about the 
formation of the Galactic disc. 

In a recent contribution, P10 showed that internal ro- 
tation in the progenitor dwarf might help to explain the 
origin of the bifurcation in the leading tail of the Sgr 
stream, as well as the puzzling shape and orientation of 
the MW dark matter halo derived by LM10. In partic- 
ular, P10 found that the same mechanism that gave rise 
to a bifurcated leading tail, i.e. a transfer of internal an- 
gular momentum into the tidal stream, could have also 
altered the orientation of the leading tail with respect 
the progenitor's orbit, directly affecting the constraints 
on the halo shape derived from stream models that as- 
sume a pressure-supported progenitor galaxy. To check 
this scenario, these authors proposed a kinematic survey 
of the Sgr remnant core, which according to their models 
should exhibit a clear signature of internal rotation. 

Here we have conducted a spectroscopic survey of the 
inner regions of the Sgr dwarf using the AAOmega in- 
strument on the AAT. Our survey covers 6 fields with a 
total area of 6ir ~ 18.8 deg 2 and provides radial veloci- 
ties for more than 1800 Sgr dwarf member stars with a 
typical accuracy of ~ 2 km/s. We have compared our 
velocity measurements against numerical N-body simu- 
lations that adopt rotating (P10) and pressure-supported 
(LM10) dwarf N-body models. We find that models with 
little or no intrinsic rotation provide a better match to 
the mean line-of-sight velocity in each of the surveyed 
fields, hence refuting the scenario proposed by P10, but 
fail to reproduce the observed line-of-sight velocity dis- 
tribution, which appears considerably more leptokurtic 
(or peaked), especially in the inner- most fields. 

Interestingly, the absence of intrinsic rotation may be 
consistent with recent theoretical models that propose 
tidal stirring as the driving mechanism for the forma- 
tion of the Sgr dwarf (Lokas et al. 2010). In this sce- 
nario discs are transformed into spheroids through a com- 
plex process that involves tidal mass loss, bar forma- 
tion/disruption and ram stripping. However, although 
these models are able to reproduce the elongated shape 
of the Sgr dwarf, as well as the absence of intrinsic angu- 
lar momentum in the remnant core, they cannot explain 
the presence of a bifurcated tidal tail with the observed 
characteristics. 

Further theoretical work is clearly necessary in order 
to reproduce the internal kinematics of the remnant core, 
as well as the locations and velocities of the associated 
stream, simultaneously. The kinematic data presented 
in this contribution will provide useful constraints for 
future numerical models that treat the dark and baryonic 
material as separate dynamical components (e.g. Lokas 
et al. 2010). Meanwhile, the origin of the prominent 
bifurcation in the distribution of tidal debris of the Sgr 
dwarf stands as an open - and enigmatic - question. 
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